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ABSTRACT 



The chariging of dust grains in astrophysical environments has been inves- 
tigated with the assumption these grains are homogeneous spheres. However, 
there is evidence which suggests many grains in astrophysical environments are 
irregularly-shaped aggregates. Recent studies have shown that aggregates ac- 
quire higher charge-to-mass ratios due to their complex structures, which in 
turn may alter their subsequent dynamics and evolution. In this paper, the 
charging of aggregates is ecamined including secondary electron emission and 
photoemission in addition to primary plasma currents. The results show that 
the equilibrium charge on aggregates can differ markedly from spherical grains 
with the same mass, but that the charge can be estimated for a given environ- 
ment based on structural characteristics of the grain. The "small particle effect" 
due to secondary electron emission is also important for determing the charge of 
micron-sized aggregates consisting of nano-sized particles. 

Subject headings: dust, extinction — interplanetary medium — ISM:general — solar 
system: general 
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Introduction 



Dust, an ubiquitous component in the universe, plays an important role in the 
thermodynamics and chemistry of the interstellar and intergalactic medium, inte rstellar 



gas dynamics, and the formation of stars, planets and planetesimals ( jJones 



19971) • Cosmic 



dust grains also interact with electroma gnetic radiation nearby, altering the observed 



spectra of remote objects (jPraine 



20031 ). The study of cosmic dust has steadily gained 
attention as technological advances makes in-situ measurements and sample-collection 
within the heliosphere feasible. Various space missions have been conducted to probe and 
investigate the composition, size distribution and structural cha racteristics of interst ellar 



and interplanetary dust, such as Ulysses, Cassini and Galileo (lAltobelli et al 



2003 



Kriiger et al. 



2010al Jbl). These dust grains provide an excellent window into interstellar 
and interplanetary processes, carry information on the origin and evolution of their parent 
bodies, and reveal the intrinsic properties of the environments where these grains originate. 
The in-situ measurements also provide opportunities to test and validate various theories 
related to cosmic dust. 

Dust grains in the solar system environment are subject to three charging processes 



(IMendis and Rosenberg 



1994 ) — for dust grains close to the Sun, strong UV radiation 
can excitep and liberate photoelectrons from the surface of the dust grain and charge 
the grain positively. At the same time, free moving electrons and ions in the solar wind 
constitute plasma charging currents , although the strength of the currents may be less than 



the photoelectric charging current (IKimura and Mann 



1998a). Between the termination 



shock and the heliopause, the boundary of the heliosphere, the plasma temperature rises 
to 2 X 10® K as a result of the constant interaction between the interstellar medium and 
the solar wind. Secondary electron emission then becomes the dominant charging process 



(IKimura and Mann 



1998ah . 
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Ascertaining the charge on cosmic dust grains in the solar system is essential for several 
reasons. First, charged grains are subject to the Lorentz force, which can significantly alter 
their trajectories. When performing on-board measurements and sample-collection, the 
grain dynamics must be taken into account to identify the origin of the measured/captured 
grains. Second, the high plasma temperature between the termination shock and the 
heliopause highly charges interstellar dust grains entering the solar system. The smaller 
mass grains are affected most strongly by the Lorentz force caused by the magnetic field 
at the heliopause. As a result, their trajectories are significantly altered and these grains 
have much smaller probability of reaching the inner solar system. Space missions devoted 



to studying t he interstellar dust in the inner so l ar system have to ta . 



into account ( Frisch et al. 



1999 



Landgral 



2000 



Linde and Gombosi 



^e th is filtering effect 



20001). Third, grain 



charging theory has successfully explained many astrophysical phenomena, such as the 



mysterious spokes of Satu rn's B ring, and the 



above the moon's surface (IMitchell et al 



2006 



^unar "horizon glow " due to dust levitating 



Stubbs et al 



20061). 



C harging of cosmic dust g r ains has been examined in detail due to the reasons listed 



above (IFeuerbacher et al. 



Kimura and Mann 



1973 



Draine and Salpeter 



1979 



Mukai 



1981 



Chow et al. 



19931 : 



1998al ). However, most of the previous work assumes a simple geometry 
for dust grains, such as a homogeneous sphere. Dust grains naturally occurring in space 
constantly go through formation processes such as nucleation, condensation, coagulation 
and destruction, and as a result they assume more complex structures such as ellipsoids 
or fluffy aggregates. Aggregates are a common structure for interplanetary, cometary and 
interstellar dust. The data and samples collected from space have shown that a porous fluffy 
structure, consisting of many small subunits, can be found among interplanetary dust grains 



and in cometary bodies (IBrownlee et al. 



1980 



Greenberg and Hage 



1990 



Hu and Winarski 



20111 ). Although no direct evidence shows that interstellar dust grains assume the same 



structure, data analysis from remote observation supports the existence of aggregate 



( Jones 


1988; 


Woo et al. 


1994; 


Wurm et al. 


2003) 



The charging of aggregate d ust grains has been studied recently both in la 



astrophysical environments ( iWiese et al. 



20101 : 



Ilgner 



2012 



Okuzumi et al. 



Doratory and 



201l[ ): these 



results have shown that aggregates tend to acquire more charge when compared to spherical 
grains of the same mass due to the porous/fluffy structure of the aggregate, and that 
charged aggregates have a significant effect on subsequent dust evolution. However, there 
has not been a detailed study of characterizing the charge on aggregates based on structural 
characteristics. 

In this paper, a 3D model is employed to calculate the charge on aggregate dust grains 
under charging conditions particular to solar system environments. Three different charging 
processes are taken into account — plasma charging currents, secondary electron emission 
and UV radiation. For charging processes dominated by secondary electron emission at the 
heliopause, we show a charge enhancement for all the aggregates compared to spheres of 
the same mass. This charge-to-mass ratio is significantly higher for aggregates consisting of 
nano-sized grains as a result of the small particle effect. For photoemission, the charge on 
aggregates may be more or less than the charge on spheres with the same mass, depending 
on the magnitude of the variable solar UV photon flux. Charge estimate models are 
proposed for both secondary electron emission and photoemission, and it is shown that 
the charge on aggregates can estimated based on structural characteristics such as the 
compactness factor. 



2. Charging Model 



The charge on a dustl grain embedded in plasma is determined by 
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where is the current contributed by the j charging process. The charge on the grain 



reaches equihbrium when = 0- isolate d spherical g rain, Equation [T] can be 



solved analytically to yield the equilibrium charge (iGoertz 



19891). However, this highly 



idealized circumstance is almost never satisfied in astrophysical environments. As discussed 
above, the cosmic dust grains often assume irregular shapes, thus require numerical 
simulation. The charging current density for three charging processes — plasma charging, 
secondary electron emission and photoelectric emission, are given below. Silicates have been 



identified as one of the major constituents for cosmic dust grains (ISavage and Mathis 



McCarthv 



1979 



19801 ). thus silicate grains with a density of 3.2 g cm ^ ( iDraine and Salpetei 



19791 ) are used as the grain material in this study to demonstrate the different charging 



processes. 



2.1. Collection of Plasma Particles 



The current density to a spherical grain can be found from Orbital Motion Limited 
theory (OML), which is based on the conserv ation of energy and angular momentum 



( IWhipple 



1981 



Laframboise and Parker 



19731 ). The current density to any point on the 



surface of a grain due to the collection of a given species of plasma particles is given by 




Vsf{vs) COS adv 



(2) 



where Ug and qs are the number density and charge of the given species, Vs is the s peed of 



the p articles, f{vs) is the distribution function which is assumed to be Maxwellian (IGoertz 



19891 ). and a is the angle between the impinging velocity and the surface normal of the 
dust grain. In the three dimeri sional case, we use spherical coordinates {v, 6, 0) in v space 
( iLaframboise and Parker! 1 19 731 ). The differential velocity dv^ can be written as 



dv^ = v^dvdVt, 



(3) 
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This allows the integration over the speed to be separated from the integral over the open 
solid angles, allowing Equation |5] to be rewritten as 



oo 

3 . 



Js = nsQs / Vsfi^s)dvs / / COS ad^. (4) 
The integration over speed is easy to carry out with Vmin given by 

0, qsf > 

(5) 



q^ip < 



where ip is the surface potential of the grain and is the mass of the plasma particle. 
For a point on the surface of an isolated sphere, the integral over the open solid angles (a 
hemisphere) is also simply evaluated. However, on the surface of an aggregate, not all of 
the incident angles are open to the incoming particle flux. Thus, the differential solid angle 
dQ requires numerical simulation for aggregates, which is discussed in Section 2.4. 



2.2. Secondary Electron Emission 



Energetic primary electrons can release secondary electrons from the surface of a grain 
upon impact, which constitutes a positive charging current. It has been shown that the 
secondary electron yield is enhanced when the dimensions of the monomers ar e comparab 



to th e primary electron penetration depth, the so-called small particle effect (I Chow et al 



19931 ). Since the size of a representative interstellar dust grain is normally less than 10 /im. 



we employ a model which takes the small-particle effect into account in determining the 



yield , S, as a function of Eq, the initial energy of the primary electron ( iDraine and Salpeter 

mm, 



8Eo/E„ 



[l + E^/E^f 



exp 



-4a 



(6) 
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Here 



fiix) 

f2{x) 



1.6 + 1.4x^ + 0.54x^ 

1 + 0.54x4 
l + 2x^ + 



(7) 



and a is the radius of the grain. The maximu m yield (5m , and the corresponding maximum 



ener gy Em., are 2.4 and 400 eV for sihcates (iMukai 



nm (jPraine and Salpeter 



198l[ ). The escape length A is 2.3 



19791). The projected range R gives the penetration depth of a 
primary electron into matter alon g the incident direction, an d is determined based on Eq 



as shown by Draine and Salpeter (jPraine and Salpeter 



19791 ). 



Thus, the current density due to secondary electron emission is calculated as 




vf{v)cosa6{Eo)dv'^x / p{E)dE, 
where p{E) is the energy distribution of the emitted electrons. It can be written as 



P{E) 



E 



1 1^ E 
1 + - 



-3/2 



(9) 



19891) 



2 \ kTgQf^ 

where Tgec is the temperature of the emitted electrons and is set to be 2 eV (iGoertz 
The lower limit of the integral is Emin = max{0, eip), with being the surface potential of 
the target grain. Equation [3] can be used in [8] to yield 



J. 



vy{v)6{Eo)dv II cos adn / p{E)dE, 

J E 



(10) 



which has a form similar to that of Equation HI with the only term dependent on the 
geometry of the aggregate being the integral over dVL. 



2.3. Photoelectric Emission 

Incoming photons with energy hv > W, the work function of the material, can excitep 
and liberate electrons from the surface, and thus constitute a positive charging current. 
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Assuming an isotropic source of UV, the photoelectric current de nsity can also be separate d 



into integration over the photon energy and the incident angles (IKimura and Mann 

POO 

Jph = Qe I Qabs{hv)F{hv)Y{hv) d{hu) 



1998af ): 



w 



cosadn I f{E)dE, (11) 



where hu is the photon energy and F{hv) is the photon flux at a given distance from the Sun, 
which can be easily obtained through satellite measurements. The absorption efficiency Qabs 
depends on the grain radius a and the wavelength of the incoming photon, A. If 2TTa > A, the 
interaction of the grain and the photon can be regarded as elastic scattering and Q nh.^ = 1- 



If 27[a < A, Mie scattering is often used to obtain Qabs 



Bohren and Huffman 



work function of silicates, W, is 8 eV based on er npirical value ( 



The photoelectric yield Y{hv) is estimated using (jPraine and Salpeter 



1983h. The 



Draine and Salpeter 



19791): 



19791 ) 



Y{hu) 



[hu -W + er 



)2 _ g2 . 



1-fl 



-f 

a 



(12) 



wher e the escape length of the photoelectrons le is 1 nm for silicates (jPraine and Salpeter 



19791 ) . emin is the minimum energy needed 



6 eV for silicates (jPraine and Salpeter 



br photoelectric emission to occur, and is set to 



19791). 



The energy distribution of the photoelectrons is given by f{E) and must be taken 
into account if the potential distribution about the surface is such that some of the 
photoelectrons return instead of escaping. Both laboratory and space experiments indicate 



that the photoelectrons are emitted isotropically with a Maxwellian distr i bution 



br the 



energy at a characteristic temperature of 1-2 eV (jHinteregger et al. 



1959 



Grard 



Wrenn and Heikkila 



19731 : 



19731 ). Thus, the energy distribution f{E) of the photoelectrons is 
E ( E 



f{E) 



exp 



(13) 



{kTpfi)'^ \ kTph 

with temperature Tph of the photoelectrons set to be 1 eV in the current study. The lower 
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limit Emin of the integration is max(0, eyjs). Thus the integration yields 



1, qey:>s < 



exp(— ^^), qe<fs>0- 



With the aid of Equation [H] the photoelectric current density can be written as 

Jph = qe Qabs{hv)F{hv)Y{hu)d{hv) / / cosarffi, geV^s < 

Jph = Qeexp^-'^^^ j Qabs{,hv)F{hv)Y{hv)d{hv) jj cosadQ qeVs>^- 



(15) 



Once again, the current density depends on the aggregate geometry through the 
integral over the open angles. Since only the side of an aggregate currently facing the sun 
is illuminated, the photon current is estimated by dividing by a factor of two (assuming 
isotropic flux). This should give an upper bound for the photoelectric current as the 
rotational period of a micron-sized dust grain is very short compared to the equilibrium 
charging time (see Section 3.2). 



2.4. Line-of-Sight Approximation 



The charging code OML_LOS calculates the electron and ion fluxes by determining 



the open lines of sight (LOS) to the 
detailed description can be found in 



points on the su r 



Matthews et al. 



ace of each constituent monomer. A 



(120121 ): here a brief summary is given. 



Electrons and ions coming from the surrounding plasma are assumed to move in straight 
lines and are captured at the points at which their straight line trajectory intersects a 
monomer, as illustrated in Figure [1] . 

The surface of each monomer is divided into many equal-area patches. Test 
directions t from the center of each patch (the so called lines of sight) are determined 



Fig. 1. — Open lines of sight to given points on a monomer in an aggregate are indicated by 
the shaded regions. Charging currents to a given point are only incident from these directions. 
The dotted line indicates an emitted electron which is recaptured by another monomer along 
a closed line of sight, while the dash-dotted line indicates an emitted electron that escapes 
along a free line of sight. 

to be blocked if they intersect any other monomer in the aggregate, or the monomer 
in question (LOS^ = 0), and open otherwise (LOS^ = 1). The integration over the 
angles in Equations HI [10] and [15] is replaced by the line-of-sight factor, which is equal 
to the sum of the open lines of sight multiplied by the cosine of the angle of the test 
direction with respect to the surface normal and by the area of the patch on a unit sphere, 
LOS = ff cos adQ = LOSi cos atA(cos 6')A0. The net current of species s to a given 
patch at a given time, Is{t), is found by multiplying the current density by the area of the 
patch. A: Is{t) = Js{t)A. Summing over species s provides the change in the surface charge 
on the patch during a time interval dt, dQ{t) = '^Is{t)dt. The contribution to the dipole 
moment is given by dD{t) = ^/s(t)R(it, where R is the displacement vector from the 
patch to the center of the grain. The current density Js{t) depends on the potential at the 
center of the given patch, which in turn depends on the charge and dipole moment on each 
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monomer, 

where Qi and Dj are the charge and dipole on the zth monomer, and Tj is the distance 
vector from the center of the i^^ monomer to the patch. The solution requires numerical 
iteration until equilibrium is reached. The change in the charge and dipole moment of 
each monomer is then obtained by adding the contribution of all the patches. The change 
in the charge and dipole moment of the aggregate is obtained by adding the contribution 
from each of the monomers. This process is iterated in time until the average change in 
aggregate charge becomes neghgible, dQagg < O.OOOlQagg, at which point the net current to 
the aggregate will be near zero. 

In computing the current due to secondary electron emission or photoemission, an 
electron escapes from the aggregate only if the randomly chosen escape direction is along an 
open line of sight. Electrons which are released along a blocked line of sight are recaptured 
by another monomer within the aggregate, leaving the total charge of the aggregate 
unchanged, but the charge distribution on the surface is altered. 



2.5. Aggregate Builder and Compactness Factor 



The numerical code Aggregate Builder was used to create aggregates through the 
coagulation of spheres using a com bination of particle-cluster aggregation (PCA), a nd 



cluster-cluster aggregation (CCA) (IMatthews et al. 



2007 



Matthews and Hyde 



20091). 



During PCA, a target particle is placed at the origin, and a single particle is released 
at the boundary of the simulation box with its velocity directed towards the center of 
target particle plus an offset. A successful collision is detected if constituents of the 
target and projectile actually touch or overlap. The grains are assumed to have relative 
velocities that are too low for any restructuring to occur, and to stick at the point of 
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contact (jWurm and Bluirull998l : iBlum and Wuriru |2000| ) . New aggregate parameters are 
then calculated, and the resultant aggregate is saved to a library. In the case of CCA, 
small aggregates from the previously saved library are employed as the target grains, with 
the incoming grain either a spherical monomer or an aggregate randomly selected from the 
same library. 

While the structure of the aggregates (characterized by the compactness factor, 
described below) depends on the plasma environment in which it grows, the charge on 
an aggregate within a given environment is in turn a function of the compactness factor 



( iMatthews et al. 



2OI2I ). Thus a large number of aggregates were built covering a wide 
range of compactness factors assuming a neutral environment. The aggregates from the 
library were then charged through OML_LOS using the parameters representing different 
astrophysical environments. 



The compactness factor, $0-, defined by Paszun and Dominik (jPaszun and Dominik 



20091 ). is used to characterize the structure or fluffiness of an aggregate consisting of 
spherical monomers. 



(17) 



where is the number of monomers in the aggregate, a is the constituent monomer radius, 
and Ra is the radius of the average projected surface area, defined as 

Rn 



with cr being the projected surface area averaged over many orientations. Figure H] shows 
a representative aggregate consisting of mono-disperse monomers, with R„ and the outer 
radius Rmax indicated. For compact aggregates, the volumes of the two spheres with these 
radii are approximately equal. For open aggregates, the ratio of the two approaches zero. 

Figure [3] shows the log-log plot for the mass of aggregates consisting of monodisperse 
monomers of different radii (5 nm < a < 500 nm), ranging in size from two to 200 
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Fig. 2. — Illustration of the compactness factor. The inner shaded area corresponds to a 
sphere with radius R^, the outer shaded area to a sphere with maximum radius, Rmax- 

monomers, as a function of the compactness factor. For each of the groups, as the mass 
increases, the compactness factor decreases, indicating a fluffier structure. This power-law 
relationship confirms that the compactness factor can be used to characterize the structure 
of the aggregates. 



3. Results 

3.1. Charging of Aggregates with Secondary Electron Emission 

In this section, the charging of aggregates is examined by including both plasma 
currents and secondary electron emission. The aggregates are charged using parameters for 
conditions in the heliosheath, the region between the termination shock and the heliopause. 
The plasma temperature in this region rises to 2 x 10^ K, a result of the constant interaction 
between the solar wind and the interstellar medium. Thus, secondary electron emission i s 



the do minant charging process due to the high electron temperature (iKimura and Mann 



1998aJ). The plasma density and temperature at 150 AU are rif. = rii = 2 x 10^ m ^ and 
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Fig. 3. — Mass of aggregates as a function of the compactness factor. The radius of the 
constituent monomers is given in the legend. 



Te = Ti = 2 X 10^ K flSchwenn 



1990 



Pauls and Zank 



19961 ). It is shown that the collective 



charge on an aggregate consisting of nano-sized grains is appreciably enhanced due to the 
small particle effect on each subunit. Two models for approximating charge on aggregates 
are proposed and the charge-to-mass ratio of the aggregates is compared to that of the 
spheres with the same mass. 



3.1.1. Charging Time 

Before estimating the equilibrium surface charge on the aggregates, the time to reach 
the equilibrium condition needs to be considered, for depending on the plasma parameters 
and the dynamic processes being considered, the equilibrium condition is not always 
satisfied for grains of all sizes. 

The dominant current determines the polarity of the equilibrium charge, while the 
non-dominant current determines the charging time, Teq. As the grain charges, the relative 
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contribution of the non- dominant current increases to balance the dominant current. Due 
to the high temperature of the plasma near the heliopause, secondary electron emission is 
the dominant charging process, determining Q. Thus be approximated by \Q/I(,\. 

Generally, Teq increases with decreasing dust radii a, approximately according to r^q oc a~^. 
Figure m shows the charging history of a dimer consisting of two 5 nm-radius monomers, the 
smallest aggregate in the simulation. The maximum charging time is approximately 1 x 10^ 
s, which is less than 5.75 x 10^ s, the time needed for interstellar dust grains to travel 1 
AU with a constant speed of 26 km/s. Since th e typical distance between the heliopause 



and termination shock is 50 AU (ISchwenn 



1990l ). all the aggregates in the simulation are 



assumed to reach equilibrium within traveling a distance of 1 AU. 




Fig. 4. — Charging curve for a dimer consisting of two 5 nm-radius monomers. Teq is 
determined by the point where the absolute change in the charge is less than 0.1% of the 
equilibrium charge. 



3.1.2. Model for Estimating Aggregate Charge 

The equilibrium surface charge on aggregates is plotted using both the number of 
monomers (Figure [5] ) and the compactness factor (Figure [6^). The aggregates in each 
group consist of up to 200 mono-disperse monomers with radii ranging from 5 nm to 500 
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Number of Monomers 



Fig. 5. — Surface charge on aggregates as a function of the number of constituent monomers. 
The hnear fits have the same slope for all monomer sizes with standard error less than 2%. 
The radius of the constituent monomers is indicated in the legend. 



nm. Figure 5 shows that the aggregate charge is related to the number of monomers by 



Qagg OC N^'^'^ 



where N is the number of monomers within an aggregate. 



(19) 



The aggregate charge can be predicted based on both the number of the monomers 
or the compactness factor. However, it is difficult, if not impossible, to determine the 
number of monomers within an aggregate measured in situ, while compactness factor can be 
obtained through remote observation. In fact, much research has bee n devoted to relatin g 



the m o rphological st r ucture of aggreg ates to their optical properties ( 



1998b: 


Kimura 


2001; 


Shen et al. 


2008) 



Kimura and Mann 



factor may serve as a useful tool when investigating the dynamics of interstellar dust grains 
in the outer heliosphere. 



The equilibrium charge as a function of the compactness factor is shown in Figure [6]^a) 
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Fig. 6. — (a) The charge on aggregates consisting of mono-disperse monomers of different 
radii, as indicated by the legend, and (b) the surface charge divided by the capacitance of 
a single monomer. The small particle effect is clearly evident for aggregates composed of 
monomers of size a = 5 nm. 

for aggregates consisting of monodisperse spheres with different radii. Each group can be 
fit with a straight line of the same slope on a log-log plot with charge related to the number 
of constituent monomers by 

Qagg (X ^-^-^ (20) 

The results clearly demonstrate that for each group, the surface charge on the aggregate 
increases as the fluffiness of the aggregate structure increases, with the surface charge of a 
sphere being the lower limit as $0- approaches one. This indicates that using the surface 
potential of a sphere of an equivalent mass to calculate the charge on an aggregate leads to 
charge underestimation. 

In Figure [6]^b), the charges on the aggregates are divided by 47reoa, the capacitance of 
a single monomer, to yield the effective potential for the aggregates. After eliminating the 
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monomer size factor, it can be seen that aggregates consisting of monomers with a > 10 nm 
fall on the same line, while the aggregates with the smallest monomers, a = 5 nm, exhibit 
a substantially higher y-intercept. The aggregates consisting of monomers with a = 10 nm 
lie in between. It is evident that the collective contribution of the higher potential achieved 
by nm-radius grains (within the much larger aggregate) caused by the small particle effect 
is significant, and needs to be taken into account when estimating the charge on aggregate 
structures. 



* 5 nm 
■ 10 nm 
o 50 nm 
+ 100 nm 
□ 500 nm 



10 10 10 

Mass (kg) 



Fig. 7. — Comparison of the surface potential on aggregates (data points) and spheres having 
the same mass (solid line). The monomer radius within the aggregates is indicated by the 
legend. 



In a given plasma environment, spherical grains with a > 10 nm reach the same 
equilibrium potential, independent of their radii, while spheres with radii a < 10 nm 
have a greater pot ential due to the small particle effect from secondary electron emissions 



(IChow et al. 



19931 ). However, as shown in Figure [7], the surface potential of an aggregate 
clearly does not follow this trend. The surface potential of an aggregate in this case is 
defined as 

Qagg 



(21) 



0' mass 



where Qagg is the total charge on the aggregate and Vmass is the radius of a solid silicate 
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sphere having the same mass as the aggregate. Overall, the surface potential of aggregates 
shows greater fluctuation and is generally greater than that of a sphere with the same 
mass, due to the greate r surface area of t he aggregate. A similar effect has been seen in 



an experimental study (IWiese et al. 



20101 ). Aggregates consisting of monomers with a = 5 



nm and 10 nm have a surface potential which is significantly higher than that for a sphere 
with the same mass. This is caused by the high positive charge each co nstituent monomer 



carries as a result of the small particle effect, consistent with prediction (IKimura and Mann 



1998af ). 



3.2. Charging of Aggregates with UV Radiation 



The typical plasma environment at 1 AU is used is used to illustrate charging 
by photoemission, as phot oemission is the dominant charging process at this distance 



(IKimura and Mann 



1998al ). Plasma and UV radiation parameters vary greatly over spatial 
distance and with time. However, the current purpose is to demonstrate the charging 
of aggregate grains compared to spherical grains, so more emphasis is placed on the 
characteristics of aggregate charging rather than modeling a specific environment. The 
pla sma density an d temperature at 1 AU are ng = = 6 x 10^ m~^ and Te = Tj = 2 x 10^ 



K (jSchwenn 



19901 ). Only electrons and singly ionized hydrogen are consi dered, with oth er 



plasma components neglected due to their relatively small contribution (jSchwenn 



199o[ ). 



Instead of simulating the photon flux in specific regions and time periods, the product of the 
yield and solar flux integrated over the spectrum is left as a free parameter and varied from 
5 X 10^^ m~^s~^ to 1.5 X 10^^ m~^s^^, within the photoemission current densities expected 



at 1 AU (IWhipple 



19811 ). The constituent monomers are taken to be silicate grains with 
radii a = 50 nm and 1 fim. The absorption efficiency Qabs is set to be unity, as the grain 
radius greatly exceeds the photon wavelength. The small particle effect for photoemission 
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is also neglected due to the large radius of the grains. 

In Figure [SI the time evolution of the aggregate surface charge is compared to that of 
an equivalent sphere for three different photoemission current densities of 6 x 10^^ m^^s~^, 
9 X 10^^ m~^s~^ and 1.2 x 10^^ m^^s~^, respectively. The results indicate that aggregates 
and spheres may have charges of opposite polarity under the same conditions (Figure [Hl^c)), 
and may be either more or less highly charged than an equivalent sphere depending on the 
magnitude of the photoemission current density compared to the plasma current density, 
(Figure IH](a) and (b)). This is due to the porous structure of the aggregate. A highly 
irregular object has a greater surface area and is thus a ble to absorb more of the emitted 



20101 ). By the same 



electrons, as also shown in a recent experimental study (IWiese et al. 
token, when the photoemission current is very strong, the porous aggregate has more 
surface exposed to the UV photons, yielding a greater positive charge. 



3.2.1. Model for Estimating Aggregate Charge 

Again, we characterize the equilibrium charge on the aggregates due to plasma and 
photoelectric charging using both the number of monomers and the compactness factor. 
Using a photoemission current density of 1.1 x 10^^ m~^s~^, the charge on aggregates is 
plotted as a function of the number of monomers and the compactness factor in Figure [91 
Based on Figure [HI (a), the charge can be estimated as a function of the number of 
monomers, 

Qagg CX N''-^\ (22) 

where N is the number of the monomers within an aggregate. Figure [9] (b) indicates 
that a linear relationship on a log-log scale can also be obtained for the charge and the 
compactness factor, 

Qagg (X (23) 
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Fig. 8. — Evolution of charge on aggregates compared to that of an equivalent sphere. The 
photoemission current density is 1.2 x 10^^ m~^s~^ in (a), 9 x 10^^ m~^s~^ in (b) and 6 x 10^^ 
m~^s~^ in (c). The aggregate charge shown is the average for six different aggregates with 
A'^ = 13 monomers. 

The same exponential factor for both monomer sizes serves as strong evidence that 
aggregate charge is a function of the aggregate structure. The charge on aggregates can 
be estimated either by the number of the constituents, or the structural characteristics 
(fluffiness of the aggregate). While determining the charge on an aggregate based on the 
number of constituent monomers seems intuitive, the information is often hard or infeasible 
to obtain. Structure characteristics, on the other hand, can be obtained through the 
scattering and absorbtion interaction between aggregates and light. The power-law relation 
between the compactness factor and the aggregate charge also provides an indirect but a 
rather accurate method of determining the morphology of interplanetary dust. The Cosmic 
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Fig. 9. — Surface charge on aggregates as a function of (a) the number of monomers and (b) 
the compactness factor. The radius of the constituent monomers is indicated in the legend. 
The hnear fits have the same slope for two sizes with standard error less than 3%. 

Dust Analyser (CP A) on the Cassini s pacecraft has successfully detected the charge on 



interplanetary dust (IKempf et al. 



20041 ). If the composition and size distribution of these 



grains is known, along with the solar wind conditions, the structure of these grains may be 
obtained based on the charge estimate models proposed above. 



4. Discussion and Conclusion 

A numerical model has been used to calculate the charge on aggregate structures in 
astrophysical environments, including primary plasma currents, secondary electron emission, 
and photoemission. It is shown that the charge on aggregates is strongly correlated to 
structural characteristics (Figure |6] and Figure [9]d, as measured by the compactness factor). 
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In general, porous aggregates, with their greater surface area, are more highly charged than 
an equivalent mass sphere, with the sign of the charge being determined by the dominant 
charging current. The substantial increase in charge-to- mass ratio for aggregates in the 
region of the heliosphere (Figure [7]) will have a significant effect on the dynamics of these 
grains, greatly influencing the mass distribution of interstellar dust grains detected within 
the solar system by instruments such as the dust detector on Ulysses. 

It is interesting to note that the relationship between charge and structural factors 
(Equations Uni [201 1^5]) . is the same for the two different charging environments. 
This is a result of the LOS factor for an aggregate being independent of the monomer 
size within the aggregate, as long as all of the spherical monomers have the same radius. 
The relationship between charge and aggregate structural characteristics for polydisperse 
monomer populations is the subject of current research. 

Finally, the relationship between the charge on an aggregate and the charge on an 
equivalent sphere can vary greatly depending on the magnitude of the non-plasma currents. 
This is seen for the aggregates charged by photoemission in Figure [Hi The values used 
for the photoemission current density in the three cases shown are all within the range 
expected for solar UV flux at lAU, which varies greatly depending on solar activity. Thus 
the charging history of aggregates in space can vary greatly over time, and is markedly 
different from the charging history of a spherical grain. Further results exploring these 
differences will be presented in an upcoming paper. 
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